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ABSTRACT 

An  earlier  method  for  the  inverse  diffraction  problem,  using  only  the 
backs cattering  cross  section  as  a  function  of  aspect  angle,  yields  the  profile 
curve  for  a  strictly  convex  body  of  revolution.  If  the  profile  curve  has 
straight  line  segments,  the  body  is  not  strictly  convex.  Nevertheless,  it  is 
shown  that  the  method,  when  applied  to  such  a  body,  yields  the  correct 
profile  except  for  a  shortening  of  the  straight  segment  by  a  factor  which  is 
at  worst  0.7. 
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SECTION  I 


INTRODUCTION 

One  of  the  simplest  approaches  to  the  inverse  diffraction  problem  is 
based  on  the  geometrical  optics  backscattering  cross  section,  cr.  For  a 
smooth  strictly  convex  body  of  revolution,  cr  is  simply  related  to  the  curva¬ 
ture  k  of  the  profile  curve.  In  References  [  1  ]  *  and  [  2  ]  it  was  shown  that 
from  a  knowledge  of  a  as  a  function  of  aspect  angle  one  could  recover  the 
profile  curve  by  integrating  a  simple  ordinary  differential  equation.  The 
cross  section  function  a  is  discussed  in  Section  II,  and,  in  Section  EH,  we 
present  a  simplified  version  of  the  earlier  theory,  leading  to  an  integral 
formula  for  the  profile  function. 

A  body  of  revolution  whose  profile  curve  contains  a  straight  line  seg¬ 
ment  (such  as  a  cone  or  a  cylinder)  is  not  strictly  convex.  Nevertheless, 
it  was  found  experimentally  in  [  3  ]  that,  in  certain  cases,  satisfactory  re¬ 
sults  could  be  obtained  for  such  targets  using  the  method  derived  in  Section  III. 
In  order  to  understand  the  domain  of  applicability  of  the  method  we  therefore 
apply  it,  in  Section  IV,  to  the  case  of  an  arbitrary  straight  line  segment. 

This  corresponds  to  a  conical  section  of  the  target.  In  so  doing  we  use  the 
physical  optics  representation  for  <r.  This  leads  to  a  triple  integral.  Two 
of  the  integrations  can  be  carried  out  approximately  by  the  two-dimensional 
stationary  phase  method,  and  the  remaining  integration  can  be  carried  out 
exactly.  We  find  that  the  resulting  profile  reproduces  the  assumed  straight 
line  segment  except  for  an  error  in  the  axial  component  of  the  segment.  The 
error  vanishes  for  a  cylindrical  section  and  is  a  maximum  for  a  full 


^Numbers  in  brackets  denote  references  cited  on  page  17. 
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(truncated)  cone.  In  that  case  we  find  that  the  cone  is  shortened  by  a  factor 
of  approximately  0.7,  a  result  which  agrees  almost  perfectly  with  experi¬ 
mental  results  reported  in  [  3  ]  . 
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SECTION  n 


physical  optics  and  geometrical  optics 

REPRESENTATIONS  FOR  THE 
BACKSCATTERING  CROSS  SECTION 


We  consider  a  perfectly  conducting  body  of  revolution  illuminated  by 
a  plane  wave  (see  Figure  1) 


ikIX 

e 


(1) 


where  E^  is  a  constant  vector,  k  =  w/c  is  the  wave  number,  and 

=  -(cos  0,  sin  6)  (2) 

is  a  unit  vector  in  the  direction  of  incidence.  The  plane  of  the  figure  is 
chosen  to  contain  the  vector  I  and  the  z-axis,  which  is  the  axis  of  rotation. 


r 


Figure  1.  Perfectly  Conducting  Body  of  Revolution 
Illuminated  by  a  Plane  Wave 
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By  means  of  physical  optics  (the  Kirchhoff  method),  the  scattered  field 
can  be  represented  as  a  surface  integral  over  the  illuminated  portion  of  the 
target.  This  integral  can  then  be  reduced  to  a  single  integral  over  the 
profile  curve  (solid  line  in  Figure  1)  by  the  stationary  phase  method.  The 
details  of  this  calculation  are  given  in  [  4  ]  .  The  backscattered  far-field  at 
the  point  X'  =  -pi  is  then  given  by  (see  Equations  (2.2)  and  (2.13)  of  [4]  ): 


where 


E 

s 


27rp 


(3) 


(4) 


Here  X  =  X(s)  =  [z(s),  r(s)]  is  the  parametric  equation  of  the  profile 
curve,  N  =  -  n  is  the  outward  unit  vector  normal  to  the  curve,  and  s  is  the 
arclength  parameter.  For  simplicity  we  have  assumed  that  only  the  upper 
part  of  the  target  (solid  line  in  Figure  1)  is  illuminated. 

We  introduce  the  function 

/m  [  1/2  _  2ikl  •  X  ... 

g (0)  =  \r  I  •  n  e  ds  .  (5) 


Then 


E 

“s 


‘  k  11/2 
7r  sin#  J 


ikp+3i7r/4 

e 

2p 


(6) 
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We  also  introduce  the  backscattering  cross  section 


a 


(?) 


If  the  profile  curve  is  strictly  convex,  the  integral  in  Equation  (5)  can 
be  evaluated  by  the  stationary  phase  method,  which  is  summarized  in  the 
Appendix.  We  introduce  the  phase  function 


4>(  s)  =  2I-X(s)  (8) 

which  appears  in  Equation  (5).  We  also  introduce  the  unit  tangent  vector 
t  =  X(s)  =  and  the  unit  normal  vector  n  =  -  N  to  the  profile  curve. 

Then 


X  =  t  =  n 


(9) 


where  k  is  the  curvature.  At  a  stationary  point  s  for  Equation  (5), 


4>(s)  =  21  - 1  =  0  . 


(10) 


This  is,  of  course,  the  specular  point,  at  which  n  =  I.  At  this  point, 


<£=2I-X  =  2Kl-n*=2  , 


(11) 


and  the  stationary  phase  formula  [Equation  (50),  with  n  =  l]  yields 


g(0) 


prr(0)~l 

[kW  J 


1/2 


2ikl*  X(0)  +  i 7r /4 


(12) 
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Here  X(0)  =  [  z (0),  r(0)]  is  the  specular  point  and  k(Q)  is  the  curvature  at 
that  point.  From  Equation  (7)  we  now  see  that 


a 


7rr(fl) 
sin  6k  (9) 


(13) 


It  is  interesting  to  note  that  the  two  principal  curvatures  of  the  surface  of 
revolution  at  the  specular  point  are  k  and  (sin0j|/r.  Thus  the  Gaussian 
curvature  is  G  =  (/<sin0)/r  and  Equation  (13)  becomes 


cr 


(14) 


This  is  a  well-known  result  which  can  be  obtained  directly  by  geometrical 
optics.  We  have  derived  it  by  means  of  the  physical  optics  representation 
because  we  will  soon  consider  the  case  in  which  the  profile  curve  has  a 
straight  line  segment.  Then  k  =  0  and  Equation  (13)  is  not  valid.  In  such 
cases  the  scattering  cross  section  is  given  by  Equations  (7)  and  (5). 
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SECTION  m 


SOLUTION  OF  THE  INVERSE  DIFFRACTION  PROBLEM 
FOR  STRICTLY  CONVEX  BODIES  OF  REVOLUTION 


At  the  stationary  point  X  (0)  we  see,  from  Equation  (10),  that 


dz  cos  0  +  dr  sin  0=0. 


(15) 


Furthermore,  by  definition 


k(9) 


d£ 
ds  ’ 


(16) 


and 


ds 

dr 


I 


2 

tan  0  = 


cos  0 


(17) 


Hence,  from  Equation  (13), 


■nr  _  7rr  ds  dr  7rr  dr  ' 

sin0K  sin0  dr  d0  sin0cos0  d0 


(18) 


It  follows  from  Equations  (15)  and  (18)  that 


dr  =  —  sin0cos0d0, 

7rr 


(19) 


and 


dz  =  - 


sin  0  d  0 


7rr 


(20) 
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These  equations  can  easily  be  integrated  to  yield 


and 


cr  (fl'Jsinfl'cosfl'd#' 


(21) 


e 

m  *  z(%)  -  7  j  7m  sin'e'Ae'  ■  (22) 

90 


Equations  (21)  and  (22)  provide  an  explicit  parametric  representation  for  the 
profile  curve  X  =  [  z (0),  r(0)]  in  terms  of  the  scattering  cross  section 
(7(6). 


\ 
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SECTION  IV 


APPLICATION  TO  A  CONICAL  SECTION 


We  now  consider  a  body  of  revolution,  a  portion  of  which  is  a  conical 
section.  The  corresponding  profile  curve  then  has  a  straight  line  segment 
as  illustrated  in  Figure  2.  On  this  segment  the  curvature  k  is  zero  and  the 
derivation  in  Section  in  does  not  apply.  Nevertheless,  we  shall  examine 
the  consequences  of  applying  Equations  (21)  and  (22)  to  this  case,  using  for 
cr(0)  the  physical  optics  scattering  cross  section.  The  first  step  is  to 
calculate  a. 

From  Figure  2  we  see  that 

n  =  (-cos  /3  ,  -sin (3  ),  t  =  (-sin/3 ,  cos/3 )  (23) 


and 


X  =  [  z(s),  r(s)]  =  /zq,  r^  +  s(-sin/3 ,  cos/3),  O^s^L  .  (24) 


r 


Figure  2.  Profile  Curve  for  Body  of  Revolution 
with  Conical  Section 
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The  coordinates  of  the  terminal  point  are  then 


(Zl’  ri)  =  ^  Z(L)’  r(L^  =  (Z0  "  Lsin^>  r0  +  Lcos/3)  ,  (25) 


and 


I  •  n  =  cos0cos/3  +  sin0sin/3  =  cos  (0  -  (3 ) 


(26) 


We  also  see  that 


I*  X  =  -  ZqCos0  -  r^sin0  -  ssin(0  -  j3 )  .  (27) 


We  now  assume  that  for  0^  ^  0  <  0^  ,  contributions  to  the  scattering  cross  sec¬ 
tion  from  portions  of  the  target  other  than  the  conical  section  can  be  neglected. 
Then,  from  Equation  (5), 


g(0)  =  cos (9  -  (3  )e 


-2ik^z^cos0  +  r^  sin6^ 


„\l/2  -2ikssin(0-/3  ) ,  ,no. 

+  scos/3)  e  v  H 'ds,  (28) 


and,  from  Equation  (7) 


ff(0)  =  1777  Igl2  =  777  cos2  (0-/3) 


sin0 


sin0 


L  L 


^  ds  ^dt^  +  scos/3^1'  2  ^r^  +  tcos/3^ 


1/2  2ik(t-s)sin  (0-/3 ) 


(29) 


0  0 
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We  now  insert  Equation  (29)  into  (21)  to  obtain 


r 


2(0)-r  2 
v  '  0 


d0'  cos 


2 


cos  01  (r0  +scos/3^1//2  +  tcos/3^ 1^2 


ik$ 


(30) 


where 


<{>  =  2(t  -  s)  sin  (0'  -  /3 )  . 


(31) 


The  integral  with  respect  to  t  and  0'  in  Equation  (30)  can  be  evaluated  by 
the  stationary  phase  method.  Thus  we  apply  the  results  of  the  Appendix  with 
n  =  2.  First  we  note  that 


4>t  =  2sin(0'  -  (3),  <pQr  =  2  (t  -  s)cos(0'  -  /3 ) 


(32) 


Thus  there  exists  one  stationary  point  at 


t  =  s,  0'  =  (3 


(33) 


At  this  point 


$  = 


^tt  ^t0T 

L*Vt  ^0'0' 


0  2 
2  0 


(34) 
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Hence  det  3>  =  -4  ,  and  since  the  two  eigenvalues  must  have  opposite  signs, 
sig  <I>  =  0.  Thus  we  obtain,  from  Equation  (50),  for  9  >  (3  , 

L 

r "(0)  -  r^2  ~  cos/3  ^  ^s(ro  +  scos^)  =  2Lcos/3  (r^  +  L/2  cos/?).  (35) 

0 


It  must  be  noted  that  the  stationary  point  in  Equation  (33)  lies  in  the  domain 
of  integration  if  and  only  if  9  >  (3  .  If  0  =  /3  ,  the  stationary  point  lies  on  the 
boundary  and  it  can  be  shown  that  the  asymptotic  value  is  one-half  of 
Equation  (35).  If  9  <  (3  ,  the  integral  is  asymptotically  zero.  Thus 


r2(0)  ~  rQ2  .  So^9<(3,  (36) 

r 2(J3)  ~  rQ"  +  Lcos/3(r0  +  j  cos/3)  =  (rQ  +  j  cos/3 ) 2  +  (jcos/?)2,  (37) 


+  Lcos 


(3  <  0<  9 


(38) 


Having  determined  r(0),  we  now  insert  Equation  (29)  into  (22).  Thus 


L  L  9 


z (9)  -  zn  - - 

v  '  0  7T 


f  ds  f  dt  |  d0’  cos"  (0'  -(3 ) 


0  0  0, 


sin0'  ^r  +  scos/3^^2  ^r^  +  tcos/Sj^"2  e  1K<^  , 

(39) 


.1/2  ik$ 
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where  <pis  again  given  by  Equation  (31).  Proceeding  exactly  as  before  we 
obtain,  for  0  >  /3  , 


z (6)  -  z 


sin/3 
0  r(/3 ) 


ds  (rQ  +  scos/3 )  = 


L  (r0  +  |co80]. 


(40) 


If  we  insert  Equation  (37)  into  (40),  the  result  is 

z  ~  zn  -  Lsin/3  *Q  (41) 


where 

Q  -  (l  ♦  c2)-172, 

Thus, 

z (0)  ~  z 

0 

z (0)  ~  z^  -  QLsin/3 
and,  from  Equations  (36)  and  (38), 
r(0)  ~  rQ  , 

r(0)  ~  rQ  +  Lcos/3  , 


Lcos/3 

2r^  +  Lcos/3 

(42) 

0Q  —  0 <P  , 

(43) 

\ 

p  <  6*  01; 

(44) 

0Q^O<P, 

(45) 

/3<  0  . 

(46) 
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From  Figure  2  we  see  that  the  coordinates  of  the  terminal  point  of 
the  line  segment  are 


(47) 


Thus  we  see  that,  as  6  passes  through  the  specular  direction,  6  =  (3 ,  the 
vector  X(0)  =  [  z (0) ,  r(0)]  jumps  from  the  initial  point  to  a  point 
[zQ  -  QLsin/I  ,  r^  +  Lcos/3]  which  would  agree  exactly  with  the  terminal  point 
if  Q  is  equal  to  one.  Thus  the  term  e  in  Equation  (42)  represents  an  error 
which  vanishes  when  (3  =  7r/2;  i.  e.,when  the  conical  segment  reduces  to  a 
cylinder.  In  any  case  we  see  that 


(48) 


The  method  we  have  used  to  evaluate  the  integrals  asymptotically 
is  of  the  simplest  kind,  and  does  not  yield  the  smooth  transition  of  the 
profile  segment  from  the  initial  to  the  terminal  point.  This  transition 
could  be  obtained  by  means  of  an  asymptotic  expansion  which  is  uniformly 
valid  for  0  in  a  neighborhood  of  0  -  (3 .  Many  types  of  integrals  have  been 
expanded  uniformly,  but  unfortunately  the  integrals  considered  here  do  not 
reduce  to  one  of  the  standard  types.  The  uniform  expansion  of  such 
integrals  is  under  investigation. 

It  is  interesting  to  note  that  for  the  case  of  a  finite  cone  (see  Figure  2) 
r^  =  0  and  e  =  1.  Then  Q  =  ^F^=  «  0.  7.  In  this  case  (in  which  the  error 
is  a  maximum)  our  method  yields  a  cone  which  is  shortened  by  the  factor 
0.7.  This  result  agrees  almost  exactly  with  experimental  results  for  a  cone 
reported  in  Reference  [3]  .  For  a  cylinder,  (3  -  tt/2,  e  =  0,  and  there  is 
no  error. 
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APPENDIX 


THE  METHOD  OF  STATIONARY  PHASE  IN  n-DIMENSIONS 


The  method  of  stationary  phase,  discussed  in  Appendix  H  of  Reference 
[  5  ] ,  yields  the  leading  term  of  the  asymptotic  expansion  for  k  —  00  of  a  func¬ 
tion  f(k)  defined  by  an  n- dimensional  integral  of  the  form 

f(k)  =  ^  g(tx . tn)elk<Mtl . tn)  dtx - dtn  =  jg(T)elk<WT)dT.  (49) 

D  D 

s  \  is  a  point  for  which  ~~  (S)  =  0 
ny  &  tp, 

for  v  =  1 . n.  Each  stationary  point  S  in  the  interior  of  the  region  D 

yields  a  contribution  to  the  leading  term  of  the  asymptotic  expansion,  given 
by 

f(k)  ~  (?)  1  [ldet<^j>  l]  V2  S(S)  ek<P(S)  +  f-  sig  (4>v.^ ,  k-00  .  (50) 


A  stationary  point  T  =  S  =  /s^,  .  . .  , 


Here  (<t>  ,)  is  the  matrix  with  elements 

vy 

_  a  2<as) 

%•  8  t  8 1. 

V  J 


(51) 


and  sig  ($,,j)  denotes  the  "signature"  of  the  matrix,  defined  by 

m 

sig(«Vj)  =  £  Sgn  ^  ’ 

K — 1 


(52) 
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In  Equation  (52),  r^,  . . .  ,  are  the  eigenvalues  of  and  sgn  =  ±  1. 

We  assume  that  the  determinant  det  is  non-zero.  The  leading  term  of 
the  asymptotic  expansion  of  Equation  (49)  is  obtained  by  summing  Equation  (50) 
over  the  stationary  points  S  in  D. 


16 


REFERENCES 


[  l]  J.  B.  Keller,  The  Inverse  Scattering  Problem  in  Geometrical  Optics 
and  the  Design  of  Reflectors,  PGAPJ7,  146,  1959. 

[2]  J.  B.  Keller  and  R.  M.  Lewis,  Inverse  Scattering  Problem  for  Bodies 
of  Revolution,  RCA  Report,  1961. 

[3]  M.R.  Weiss,  Inverse  Scattering- Phaseless  Methods,  MTR-387, 

The  MITRE  Corporation,  March  1967. 

[4]  J.  B.  Keller  and  R.  M.  Lewis,  The  Fourier  Transform  Method  for 
the  Inverse  Scattering  Problem,  RCA  Report,  1964. 

[5]  R.  M.  Lewis,  Asymptotic  Methods  for  the  Solution  of  Dispersive 
Hyperbolic  Equations,  Asymptotic  Solutions  of  Differential  Equations 
and  their  Applications,  ed.  by  C.  W.  Wilcox,  Wiley,  N.  Y. ,  1964. 


17 


DOCUMENT  CONTROL  DATA  ■  R  &  D 


(Securify  classification  of  title,  body  of  abstract  and  indexing  annotation  must  be  entered  when  the  overall  report  Is  c  insal(ied) 


1  ORIGINATING  A  c  T »  v  1  T  Y  ( Corpora  te  author) 

The  MITRE  Corporation 

Bedford,  Massachusetts 

2m.  REPORT  SECURITY  CLASSIFICATION 

Unclassified 

2b.  GROUP 

3  REPORT  T 1  TLE 

PHASE-INDEPENDENT  METHODS  FOR  THE  INVERSE  DIFFRACTION  PROBLEM 

4  DESCRIPTIVE  NOTES  (Type  of  report  and  inclusive  dates) 

N/A 

$  *u  THORISI  ( First  name ,  middle  initiai ,  last  name) 

Lewis,  Robert  M. 

6  REPORT  DATE 

November  1967 

7«.  TOTAL  NO  OF  PAGES  7b.  NO.  OF  REFS 

21  5 

8a.  CONTRACT  OR  GRANT  NO 

AF19(628)-5165 

b  PROJEC  T  NO. 

4966 

c. 

d. 

9a.  ORIGINATOR*  S  REPORT  NUMBER(S) 

ESD-TR- 67-381 

9b.  OTHER  REPORT  NO(S|  (Any  other  numbers  that  may  be  as  si gne  d 
this  report) 

MTR-418 

10  DISTRIBUTION  STATEMENT 

This  document  has  been  approved  for  public  release  and  sale;  its  distribution  is 
unlimited. 

11.  SUPPLEMENTARY  notes 

N/A 

12  SPONSORING  MILITARY  ACTIVITY  DCpUty  IOT  SliTVeii” 

lance  and  Control  Systems,  Space  Defense 
Systems  Program  Office;  Electronic  Systems 
Division,  L.  G.  Hanscom  Field,  Bedford, 
Massachusetts. 

1  3  .  ABSTRACT 


An  earlier  method  for  the  inverse  diffraction  problem,  using  only  the  backscattering 
cross  section  as  a  function  of  aspect  angle,  yields  the  profile  curve  for  a  strictly  convex 
body  of  revolution.  If  the  profile  curve  has  straight  line  segments,  the  body  is  not  strictly 
convex.  Nevertheless,  it  is  shown  that  the  method,  when  applied  to  such  a  body,  yields 
the  correct  profile  except  for  a  shortening  of  the  straight  segment  by  a  factor  which  is 
at  worst  0.7. 


DD  ,F.°“  .1473 


Security  Classification 


Security  Classification 


1  4 

KEY  WORDS 

LINK  A 

LINK  B 

LINK  C 

ROLE 

W  T 

ROLE 

W  T 

ROLE 

W  T 

INVERSE  DIFFRACTION  THEORY 

Phase-Independent  Methods  v 

Aspect  Angle  Functions 

GEOMETRIC  OPTICS 

Inverse  Diffraction  Theory 

Profile  Function 

Security  Classification 


